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^ ' Abstract 

<N : 

The unusal hard x-ray burster GRO J 1744-28 recently discovered by the Compton 
■ Gamma-ray Observatory (GRO) can be modeled as a strange star with a dipolar magnetic 
O ■ field < 10 11 Gauss. When the accreted mass of the star exceeds some critical mass, its 

\o ■ 

crust may break, resulting in conversion of the accreted matter into strange matter and 

00 : 

| i release of energy - Sub — a flreM1 may form and expmd relati ™ tic 

The expanding fireball may interact with the surrounding interstellar medium, causing its 

kinetic energy to be radiated in shock waves, producing a burst of x-ray radiation. The 

burst energy, duration, interval and spectrum derived from such a model are consistent 

> ! with the observations of GRO J1744-28. 
• i— i . 

^ : PACS numbers: 97.80. Jp; 13.38Mh; 95.30Gv; 97.60Gb 

GRO Jl 744-28 is a new type of x-ray transient source, which was discovered on 2 
December 1995 by the Burst And Transient Source Experiment (BATSE) on the Compton 
Gamma-Ray Observatory (1). The bursts were detected up to energies of ~ 75keV with 
intervals between bursts of about 200 s initially. After 2 days the burst rate dropped to 
about one per hour (2). However, by 15 January 1996 the burst rate had increased to ~ 40 
per day. The burst durations lasted ~ 10 s. The burst fluences (25-60 ke V) range between 
1.7 x 10~ 7 and 6.8 x 10" 7 ergcm~ 2 ; the average fluence S = (2.7 ± 0.9) x 10~ 7 ergcm~ 2 . 
The position of the source is at the Galactic Center. For a distance of ~ 7.5 kpc, the 
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average peak luminosity was ~ 2 x 10 38 ergs _1 . Analysis of the BATSE data indicated 
that the source is a binary pulsar with a pulsation period of 0.467 s, a companion of mass 
between 0.22-1.0 M , and a binary orbital period of 11.8 days (3). Because the x-ray mass 
function is small (~ 1.31 x 1O _4 M ), the system must be nearly face-on to an observer from 
Earth with an inclination angle of ~ 18° (3-5). Furthermore, in order for the measured 
rotation period derivative to be consistent with the Standard accretion torque theory (6), 
the persistent luminosity of the source at its peak should be close to the Eddington limit 
(4) and the surface dipole magnetic field of the pulsar should be < 10 11 G (3-5). /,From 
the observed pulsed fraction and pulsar's x-ray spectrum, the strength of the local surface 
magnetic field is estimated to be several 10 12 G(^). In addition, the proportional counter 
array (PCA) experiment (2-60 keV) on the Rossi X-ray Timing Explorer (RXTE) (T) 
had detected GRO J1744-28 in the period 18 January -10 May 1996. The observations 
showed that following the earlier large bursts the flux dipped below the preburst level 
by up to 25%-30% and then made a slow quasi-exponential recovery back toward the 
preburst level. The observed recovery period lasted up to ~ 1000 s for some bursts, but 
most bursts recovered in a few hundred seconds. 

The properties of the hard x-ray bursts (HXRBs) from GRO J 1744-28 differ from those 
of other known high-energy bursts — x-ray bursts, soft 7-ray bursts and 7-ray bursts. 
First, the HXRBs are probably not type-I x-ray bursts (8). Thus, thermonuclear flashes in 
matter accreted onto the surface of a neutron star may not produce HXRBs. Second, the 
durations of the HXRBs are several hundred longer times than those of the three soft 7-ray 
repeaters even though these two kinds have similar repeat times and spectra. Third, the 
HXRBs are different from 7-ray bursts, because 7-ray bursts do not have fast repeat times 
and their spectra are much harder. On the other hand, the repeat times and spectra of 
the HXRBs are somewhat similar to those of type-II x-ray bursts from the Rapid Burster 
(2,8). This suggests that some accretion instability may be a mechanism for producing 
HXRBs. Very recently Cannizzo (9) studied the global, time-dependent evolution of the 
Lightman-Eardley instability, which might account for some observational features of the 
HXRBs. Here we propose an alternative model, in which a strange star accretes matter 
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from its low-mass companion. 

Strange matter, which is bulk quark matter, was conjectured to be more stable than 
hadronic matter (10). Studies showed that the existence of strange matter is allowable 
within uncertainties inherent in a strong-interaction calculation (11). Thus, strange stars 
may exist in the Universe. Strange stars only have crusts with masses ~ 10~ 5 M Q (12). 
However, the postglitch behavior of pulsars can be described by the neutron-superfluid 
vortex creep theory (13) which requires the crustal mass > 10~ 3 M Q . The conversion 
of a neutron star to a strange star may require the formation of a strange-matter seed, 
which is produced through the deconfinement of neutron matter at a density ~ 7-9po 
(where p is the nuclear matter density), much larger than the central density of a 1.4M 
neutron star with a moderately stiff to stiff equation of state (14). These two features 
suggest that strange stars may be formed in low-mass x-ray binaries (15) because when 
the neutron star in a low-mass x-ray binary accretes sufficient mass (perhaps > O.4M ) its 
central density can reach the deconfinement density and subsequently the whole star will 
undergo a phase transition to become a strange star. The phase transition from nuclear 
matter to strange matter under the condition of conserved charge rather than constant 
pressure may occur at a density as low as 2-3p (16). If so, strange stars may be formed 
during the evolution of protoneutron stars (17). Here we would remind the reader of some 
arguments against the existence of strange stars. For example, the disruption of a single 
strange star may contaminate the entire galaxy and essentially all neutron stars may be 
strange stars (18). In view of these uncertainties, we should only regard strange stars as 
possible stellar objects. 

According to the Standard accretion torque theory (ft)and Daumerie et al. (4) the 
binary system of GRO J1744-28 is nearing the end of the mass transfer phase. If this is 
true then the companion has transferred mass of > O.4M to the pulsar and now has a 
mass ~ O.22-O.5M . From the scenario proposed in (15) and the above discussion, we 
suggest that the pulsar in GRO J1744-28 is a strange star. 

We next discuss the burst mechanism. We assume that far from the proposed strange 
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star the magnetic field is purely dipolar, and the accretion flow is spherically symmetric. 
We consider a simple case: M pu i sar = 1.8M Q and -R pu isar = 10 6 cm. The Alfven radius 
for spherical accretion, obtained by balancing accretion and magnetic pressure, is given 
by R A = 1.8 x WL~ 2I1 BT cm, where L* is the total accretion luminosity of ~ 2 x 
10 38 ergss _1 (the Eddington limit), and B* is the surface dipolar magnetic field strength 
in units of 2 x 10 10 G which is close to the B* derived in (5). For an assumed dipolar 
magnetic field geometry, sin 2 m /r = constant, where m is the magnetic colatitude. Thus, 
at the stellar surface near a pole, the cross-sectional area of the accretion column is about 
A v « 1.8 x 10 n Ll /7 B^ /7 cm , and the corresponding radius of the accretion column is 
r p w 2.4 x 10 5 Ll /7 B* 2/7 cm. 

As the strange star accretes matter from its companion, strong pressure is formed at 
the base of the accreted matter near the polar cap, due to the gravitational attraction of 
the compact strange star. When this pressure exceeds the critical stress of the star's thin 
crust, the crust may break. Therefore, the condition under which a crust-breaking event 
takes place should be given by phg = p9, where p is the density at the base of the accreted 
matter, h is the height of the accretion column, g is the surface gravity, p is the shear 
modulus, and is the shear angle of the crust. should be about 1CT 3 for neutron stars 
in low-mass x-ray binaries to explain the bimodal magnetic field distribution of binary 
pulsars(15). We expect that for strange stars in low-mass x-ray binaries is close to this 
value because the stellar crust in both cases are replaced by accretion material. From 
(20), p ~ 2.5 x 10 27 dyncm" 2 , so we obtain the column density of the accreted matter, 
o — ph ~ 1.0 x 10 10 #_3 g cm -2 , where 6L 3 = 0/1CT 3 . The interval between crust-breaking 
events can be written as 

Tl « ^ ~ 2.2 x itfe^L^B: 4 ' 7 s . (1) 
M 

As the luminosity due to accretion decreases or the magnetic field decays, T\ increases. 
This time scale is consistent with the typical, observed intervals between the HXRBs. 
However, the accreted matter may diffuse away from the polar cap area A p before it 
builds up enough pressure to break the crust. We now estimate the diffusion time of the 
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accreted matter across the local magnetic field of the polar cap B s , which is assumed 
to be dominated by the multipole field component and in general much stronger than 
that of the dipolar field component. The transverse velocity resulting from collisional 
diffusion in the presence of a pressure gradient (VP) is approximated by ^ ~ 1.3 x 
10 2 ZT 8 3/2 B; 2 VP cm s 1 (21) where Z = 1 for hydrogen, T 8 is the matter temperature in 
units of 10 8 K and B s is the local surface magnetic field. In the present case, VP ~ pO/r p . 
So the diffusion velocity v d ~ 2 x 10~ 4 T 8 3/2 L~ 1/7 B 2/7 (B S /10 12 G)- 2 cius" 1 . Thus, the 
time scale for the accreted matter to diffuse over the length ~ 1 km is at least 10 9 s, which 
is much longer than T\. 

After the crust has been broken, the accreted matter will fail into the strange-matter 
core in a time ~ 1 ms. Subsequently releasing two kinds of energy: (i) gravitational 
energy ~ 2 MeV per nucleon, which is due to the movement of the accreted matter from 
the stellar surface to the base of the crust and (ii) deconfinement energy ~ 30 MeV per 
nucleon, due to the conversion of the accreted matter to more stable strange matter (22). 
The total released energy is E tot ~ 5.5 x 10 40 6^ 3 Ll^ 7 B* 4 ^ 7 ergs. Because the total volume 
of the strange-matter blobs formed during the conversion of accreted matter is rather 
small (~ 10 6 cm 3 ), most of E tot will be radiated through photons from the surfaces of the 
blobs. However, a fraction of the radiation energy may be absorbed and then reradiated 
as neutrinos which pass almost freely through the crust. Thus, it is expected that about 
half of the total released energy may be finally radiated in photons which form a fireball 
of volume ~ A p l (where l is close to the crustal thickness ~ 10 4 cm). Assuming that £ is 
the ratio of the fireball energy to the total released energy, we obtain the fireball energy 

P 7 ~ 2.8 x l0 4 % /2 9- 3 L 2 J 7 B; 4 / 7 ergs , (2) 

where £i/ 2 = 2£. Let T be the initial temperature of the fireball. By using the expres- 
sion ^aT^Apl ~ E 1 , we have T ~ 5.2 x lO 9 ^^ 1 / 4 K. The fireball must be contam- 
inated by baryons and we can estimate the amount of contamination. The radiation- 
dominated outflow begins when the radiation energy density w 7 = ^faT^ is equal to 
the gravitational energy density u g = GMp/r, or p = 8.7 x 10 5 (r /10 10 K) 4 gcm -3 . 
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/From (23), the column density for the radiation-dominated surface layer is given by 
a' « 3 x 10 8 / u f T 1 / 3 (p/10 6 gcm~ 3 ) 4 / 3 g cm -2 , where \i e is the mean molecular weight per 
electron (/i e = 1 for hydrogen). Therefore, the amount of the baryons loaded with the 
fireball is approximated by AM w cr'Ap ~1.4x 10 18 L 2 J 7 B~ 4/7 g. Thus, the ratio 

of the initial fireball energy to the rest energy of the contaminating baryons is defined as 
r] = ~ 2l£ 1 ^ 3 6_l^ 3 . The fireball will expand outward because of the large optical 
depth of photon-electron scattering. Because AM/M Q > 1.7 x 1(T 16 (£ 7 /10 41 ergs)(£^), 
most of the initial fireball energy will be converted into the bulk kinetic energy of the 
baryons during the expansion. When the optical depth becomes one, therefore, the rest 
radiation energy becomes negligibly small. Fortunately, as suggested by Meszaros and 
Rees (25), the expanding shell (having a relativistic factor r ~ rj) will interact with the 
surrounding medium and its kinetic energy will be converted into the random energy of 
the shell by shock waves and finally radiated through nonthermal processes in these shock 
waves. The timescale for radiation is approximated by (25) 

r 2 « 0.1<^r/10 2 )- 8 / V /3 s ~ mf/^Ll^B^W'" s , (3) 

where n is the interstellar density (~ lcur 3 ). This time scale is consistent with the 
typical, observed durations of the HXRBs. 

It is widely believed that electrons can be accelerated by shock waves to very high 
energy with the minimum Lorentz factor 7 min ~ (m p /m e )Y, assuming that all particles 
behind the shock waves have the same energy. If the shock waves can produce approximate 
equipartition between the magnetic field energy density and the particle energy density, 

1 /2 

then the strength of magnetic field is about B ~ 0.3rn G(25). The typical synchrotron 
photon energy emitted by electrons with the Lorentz factor 7 m i n at the observer's frame 
is about 

e p ^li^e-_r< /2 KeV. (4) 

It is usually expected that diffusion shock wave acceleration can produce a power-law 
electron spectrum, dN e /d^ ~ r y~ p , 7 min < 7 < 7 ma x, where the spectral indexp is typically 
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between 2 and 3 (26). The synchrotron radiation for electrons with such a distribution 
has a spectrum for photon energies larger than e p with a power-law form with the photon 
index of a = — (p + l)/2(27). So the theoretical value of a may be from about —1.5 to 
—2.0. This result is consistent with the observations by OSSE, a = —(2.0 ±0.6) (28), and 
the observations by RXTE, a ~ —1.3 (7). The bursting spectrum for our radiation model 
is similar to the spectrum of the persistent emission observed in (7,28). Theoretically 
the radiation processes at the surface near a magnetic pole of the accreting strange star 
are complicated, and to study these processes is beyond the scope of our paper. Some 
models studying the radiation processes of an accreting neutron star with a rather strong 
magnetic field may give power-law spectra with an index near —1.5(29). In addition, it 
should also be noted that the spectrum does not evolve during the bursting period because 
the spectral index for the electron distribution behind the shock wave during the shell's 
expansion is unlikely to be changed (30). 

We can estimate the recovery time scale as follows. When a burst occurs, the huge 
radiation pressure will push the accretion matter outward over a distance Ar ~ £x5a£l ? 
where L is the accretion luminosity and v r is the radial velocity of the accreted matter. 
After a burst, the accretion matter will fail back toward the strange star over a time 

r 3 ~-~2.4x s . (5) 

v r 

which is consistent with the typical, observed recovery time scale(7). 

We can compare other characteristics determined from eqs. (1), (2) and (3) with the 
observations. First, the RXTE observations on GRO J1744-28 between 29 January and 
26 April, 1996 indicate (7) that the data of the nonbursting flux from this source can be 
approximately fitted with the straight line flux (mcrab) = 2703.3 — 23. 0D, where D is 
the day number in 1996. If this expression can be extrapolated to December 1995, the 
ratio of the persistent flux on 5 December 1995 to the persistent flux on 30 January 1996 
is about 1.6. Because the interval time scale is proportional to L* 5 ^ 7 (eq.l), the ratio of 
the typical interval time scales for the HXRBs on 30 January 1996 to that on 5 December 
1995 is about 1.4. The observations from RXTE on 30 January 1996 (7) and BATSE on 
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5 December 1995 (2) showed this ratio is ~ 1.38. Therefore, eq.(l) is consistent with the 
observations. Second, for our model, the bursting flux is obtained by dividing eq. (2) by 
eq. (3), and this flux is proportional to L^ 21 . This means that the bursting flux is weakly 
dependent upon the persistent fhrx in agreement with the observations from BATSE (2,8) 
and OSSE (28). 

The surface radiation in the crust-breaking region during the bursts should show pul- 
sations amplitude close to that of pulsations during the nonbursting periods. This model 
implication agrees with the observations from OSSE (28) and RXTE (31). However, the 
reults of RXTE indicates that the bursting flux seems approximately linearly proportional 
to the persistent flux(5 > ^). The more detail discussion on the discrepancies between exper- 
imental results and model predictions is presented elsewhere (<?<?). According to accretion 
instability models, on the other hand, during the accretion instability, a great deal of 
matter falls onto the surface near a magnetic pole of the pulsar, and subsequently a great 
number of gravitational energies are released and an HXRB is thus formed. Following 
this picture, one should detect pulsations during the bursts amplitude of which is much 
larger than that of pulsations during the nonbursting periods. 

Finally, we want to remark that a similar strange star model is recently proposed 
to explain the soft 7-ray repeaters(5^). The key differences between these two models 
are (i)the crust cracking of the soft 7-ray repeater model results from the spin-down of 
the strange star instead of accretion, (ii)the amount of energy released from these two 
mechanisms differ by two order of magnitude and (iii)the strength of the dipolar magnetic 
field in these two kind of sources also differ by two order of magnitude. (iv)The time scales 
of energy release are difference by at least one order of magnitude. These differences make 
the magnetic energy density of the soft 7-ray repeaters stronger than the radiation energy 
density therefore the fireball mechanism cannot be developed. 
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